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Abstract A study of cicumnuclear star-forming re- 
gions (CNSFRs) in several early type spirals has been 
made in order to investigate their main properties: stel- 
lar and gas kinematics, dynamical masses, ionising stel- 
lar masses, chemical abundances and other properties 
of the ionised gas. Both high resolution (R^20000) and 
moderate resolution (R ~ 5000) have been used. 

In some cases these regions, about 100 to 150 pc in 
size, are seen to be composed of several individual star 
clusters with sizes between 1.5 and 4.9 pc estimated 
from Hubble Space Telescope (HST) images. Stellar 
and gas velocity dispersions are found to differ by about 
20 to 30 km/s with the H/3 emission lines being narrower 
than both the stellar lines and the [OIII] A 5007 Alines. 
The twice ionized oxygen, on the other hand, shows ve- 
locity dispersions comparable to those shown by stars. 
We have applied the virial theorem to estimate dy- 
namical masses of the clusters, assuming that systems 
are gravitationally bounded and spherically symmetric, 
and using previously measured sizes. The measured 
values of the stellar velocity dispersions yield dynami- 
cal masses of the order of 10^ to 10^ solar masses for 
the whole CNSFRs. 

We obtain oxygen abundances which are compara- 
ble to those found in high metallicity disc HII regions 
from direct measurements of electron temperatures and 
consistent with solar values within the errors. The re- 
gion with the highest oxygen abundance is R3-I-R4 in 
NGC 3504, 12-flog(0/H) = 8.85, about 1.5 solar if 
th e solar oxygen abunda nce is set at the value derived 
bv lAsplund et al. I (|2005h . 12-t-log(0/H)o = 8.66±0.05. 
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The derived N/0 ratios are in average larger than those 
found in high metallicity disc HII regions and they do 
not seem to follow the trend of N/0 vs 0/H which 
marks the secondary behaviour of nitrogen. On the 
other hand, the S/0 ratios span a very narrow range 
between 0.6 and 0.8 of the solar value. As compared to 
high metallicity disc HII regions, CNSFR show values of 
the O23 and the N2 parameters whose distributions are 
shifted to lower and higher values respectively, hence, 
even though their derived oxygen and sulphur abun- 
dances are similar, higher values would in principle be 
obtained for the CNSFR if pure empirical methods were 
used to estimate abundances. CNSFR also show lower 
ionisation parameters than their disc counterparts, as 
derived from the [SII]/[SIII]. Their ionisation structure 
also seems to be different with CNSFR showing radia- 
tion field properties more similar to HII galaxies than 
to disc high metallicity HII regions. 



1 Introduction 

The inner (^^ IKpc) parts of some spiral galaxies show 
star formation complexes frequently arranged in an an- 
nular pattern around their nuclei. These complexes are 
sometimes called " hot spots" and we will refer to these 
as circumnuclear starforming regions (CNSFRs). Their 
size s go from a fe w tens to a few hundreds of pc (e.g. 
iDiaz et al. 2000l ) and they seem to be made of an en- 
samble of HII regions ionised by luminous compact stel- 
lar clusters whose sizes, as measured from high spatial 
resolution HST images, are seen to be of only a few pc. 

The luminosities of CNSFRs are rather large with 
absolute visual magnitudes (M„) between -12 and 
-17 and Ha luminosities which are comparable to 
those shown by 30 Dor, the largest HII region in the 
LMC, and overlap with those shown by HII galaxies 
( Melnick et al. lll988t[Dfaz et al. II2OOGI : IHovos h Diaz " 



2 



2nnd and references therein). In the ultraviolet (UV), 
massive stars dominate the observed circumnuclear 
emission even in the presence of an active nucleu s 



Gonzalez-Delgado et al.1 (|l998l ): IColina et al. I (|2002l) . 



In many cases, CNSFR show emission line spectra 
similar to those of disc HII regions. However, they 
show a higher continuum from background stellar pop- 
ulations as expected from their circumnuclear location. 
The analysis of these spectra gives us the oportunity 
to measure the gas abundances close to galactic nuclei 
which, in the case of early type spirals, are expected to 
be amongst the highest metallicity regions. 

The importance of an accurate determination of the 
abundances of high metallicity HII regions cannot be 
overestimated since they constitute most of the HII re- 
gions in early spiral galaxies (Sa to Sbc) and th e inner 
regions of most late type ones (Sc to Sd) iDiaz I |l989 1 
without which our description of the metallicity distri- 
bution in galaxies cannot be complete. In particular, 
the effects of the choice of different calibrations on the 
derivation of abundance gradients can be very impor- 
tant since any abundance profile fit will be strongly bi- 
ased towards data points at the ends of the distribution. 
It should be kept in mind that abundance gradients are 
widely used to constrain chemical evolution models, his- 
tories of star formation over galactic discs or galaxy 
formation scenarios. Also, the question of how high is 
the highest oxygen abundance in the gaseous phase of 
galaxies is still standing and extrapolation of known ra- 
dial abundance gradients would point to CNSFR as the 
most probable sites for these high metallicities. 

Accurate measures of elemental abundances of high 
metallicity regions are also crucial to obtain reliable 
calibrations of empirical abundance estimators, widely 
used but poorly constrained, whose choice can severely 
bias results obtained for quantities of the highest rel- 
evance for the study of galactic evolution like the 
luminosity-metallicity (L-Z) relation for galaxies. 



set, also obtained with the same instrumentation had 
a lower resolution of 3.4 and 1.7 A in the blue and red 
spectral regions respectively and a wider total coverage 
from 3600 to 9650 A. These data were used to derive the 
physical conditions and the abundances of the emitting 
gas. 

Several CNSFR were observed in the spiral galaxies: 
NGC 2903, NGC 3351 and NGC 3504. Archive images 
obtained with the Wide Field and Planetary Camera 
2 and the NICMOS Camera 3 both on-board the HST 
were also downloaded in order to complement the study. 
Figure [1] shows the optical image of NGC 3351 with 
the slit positions used for the two sets of observations. 
Figure [2] shows the spatial distribution of the H/3 and 
continuum emission along one of the slit positions ob- 
served. As can be seen from these profiles, some clusters 
are dominated by continuum emission while in other 
cases gaseous line emission is clearly important. 



3 Summary of kinematical results 

Gas velocity dispersions were measured by perform- 
ing Gaussian fits to the H/3 A 4861 A and [OIII] A 
5007 A lines on the high dispersion spectra (Figure 
[2]). Stellar velocity dispersions dispersions were mea- 
sured using the CaT hues at AA 8494, 8542, 8662 Ausing 
the cross-correlation t echnique described in detail by 



Tonry &: Davis I (|l979f ) . Late type giant and supergiant 



stars that have strong CaT absorption lines were used 
as stellar velocity templates. 

For the 5 CNSFR observed in NGC 3351, stellar 
velocity dispersions are found to be between 39 and 
67kms~^, about 20kms^^ larger than those measured 
for the gas, if a single Gaussian fit is used. However, 
the best fits obtained involved two different components 
for the gas: a "broad component" with a velocity dis- 
persion similar to that measured for the stars, and a 



2 Observations 

As part of the programme to study the properties of 
CNSFR we obtained two sets of data with the Inter- 
mediate dispersion Spectrograph and Imaging System 
(ISIS) attached to the 4.2 m William Herschel Telescope 
(WHT). The first set of observations consisted of high 
resolution blue and far-red long-slit spectra covering the 
emission lines of H/3 and [OIII] and the Call triplet lines 
respectively. At the attained resolution of 0.4 and 0.7 
A respectively, these data allowed the measurement of 
radial velocities and velocity dispersions of the gas and 
stars in the regions down to 25 km s^^. The second 




Fig. 1 F606W (wide V) image centred on NGC 3351 ob- 
tained with the WFPC2 camera of the Hubble Space Tele- 
scope showing the WHT-ISIS slit positions overplotted for 
the kinematics set observations (left) and the abundance set 
observations (right). North is up and East to the left 
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Fig. 2 Spatial light profiles along one of the slit position 
in NGC 3351. The thin solid (red) line shows the total Ha 
emission, the thin broken (blue) line shows the continuum 
adjacent to Hq and the thick solid line (black) shown the 
pure Ha 



"narrow component" with a dispersion lower than the 
stellar one by about 30kms^^. These two components 
were found both in the hydrogen recombination lines 
(Balmer and Paschen) and in the [OIII] A 5007 A line. 
The narrow component is dominant in the H recombi- 
nation lines, while the broad component dominates the 
[OIII] one. Figure m shows this effect. 

CNSFR are seen to consist of several individual 
star clusters, although some of them seem to have an 
only knot, at the HST resolution. The derived masses 
for the individual clusters as derived using the sizes 
measured on the HST images are between 1.8 and 
8.7 X 10^ M0. These values are between 5.5 and 26 
times the mass de r ived f or the SSC A in NGC 1569 by 
Ho fc Filippenko ( 19961 ) and larger than other kine- 
matically derived SSC masses. Values for the dynam- 
ical masses of the CNSFRs are in the range between 
4.9 X 10'^ and 4.3 x lO^'Mg. Masses derived fr om the 
H/3 velocity dispersion under the assumption of a single 
component for the gas would have been underestimated 
by factors between approximately 2 to 4. 

The masses of the ionising stellar clusters of the 
CNSFRs have been derived from their Ha luminosities 
under the assumption that the regions are ionisation 
bound and without taking into account any photon ab- 
sorption by dust. Their values are between 8-0 x 10^ 
and 2.5 x 10^ Mq. Therefore, the ratio of the ionising 
stellar population to the total dynamical mass is be- 
tween 0.02 and 0.16. 

The SSC in the observed CNSFRs seem to contain 
composite stellar populations. Although the youngest 
one dominates the UV light and is responsible for the 
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Fig. 3 Blue (upper panel) and red (lower panel) rest frame 
normalised spectra of R6. Notice the absence of conspicuous 
emission lines in the red spectral range for this region. 



gas ionisation, it constitutes only about 10% of the to- 
tal. This can explain the low equivalent widths of emis- 
sion lines measured in these regions. This may well ap- 
ply to the case of other SSC and therefore conclusions 
drawn from fits of SSP (single stellar population) mod 
els should be taken with caution (e.g. lMc Cradv et al. 
2003). Also the composite nature of the CNSFRs means 
that star formation in the rings is a process that has 
taken place over time periods much longer than those 



implied by the properties of the ionised gas. 



4 Abundance analysis 

For 12 CNSFRs in the galaxies NGC 2903, NGC 3351 
and NGC 3504 chemical abundances of O, N and S 
were derived from the lower resolution spectra. All the 
spectra are characterized by very weak [OIII] AA 4959, 
5007 A Hues. Fi gure [5] shows the spectrum of a typical 
region. 

This low excitation is typical of high abundance re- 
gions in which the cooling is dominated by [OH] and 
[OIII] emission lines. In these cases the detection of 
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Fig. 4 Sections of the normalised spectrum of region R2 in 
NGC 3351. Left: H/3 and right [Oiii] A5007A . The dashed- 
dotted hne is the broad component, the dotted line is the 
narrow component and the dashed line is the sum of both 

the auroral [OIII] A 4363 A line which woukd allow 
the derivation of the electron temperature is virtually 
impossible and empirical methods based on the calibra- 
tion of strong emission lines have to be used. The most 
commonly used abundance indicator, O23, involves the 
sum of the [Oil] and [OIII] emission lines, supposed to 
be almost independent of geometrical effects in the neb- 
ulae. However, this indicator is not well calibrated at 
the high metallicity end where the oxygen lines are so 
weak that the value of the R23 practically saturates. 

The sulphur lines, much stronger than the oxygen 
lines in high metallicity regions, provide an alterna- 
tive way to derive their abundances. We have devel- 
oped a new method for the derivation of sulphur abun- 
dances based on the calibration of the [SIII] electron 
temperature vs the empirical parameter SO23 defined 
as the quotient of the oxygen and sulphur abundance 
parameters O23 and S23 and the further assumption 
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Fig. 5 Blue (upper panel) and red (lower panel) spectra 
of region R12 in NGC 2903 



that T([SIII]) ~ T([SII]). Then the oxygen abundances 
and the N/0 and S/0 ratios can also be derived. Figure 
[5] shows this calibr ation for wh i ch da ta from different 
sources as listed in lPiaz et al. I ( 2007 ) have been used. 
The actual fit to the data gives: 



te{[SIII]) = 0.596 - 0.283logSO2 



0.199{logSO23)'^ 



Only for one of the regions, the [SIII] A 6312 A line 
was detected providing, together with the nebular [SIII] 
lines at AA 9069, 9532 A , a value of the electron temper- 
ature of Te([SIII])= SAOOttllo^- This value is slightly 
higher than predicted by the proposed fit and is repre- 
sented as a solid upside down triangle in Figure [H] 

The derived oxygen abundances are comparable to 
those found in high metallicity disc HII regions from 
direct measurements of electron temperatures and are 
consistent with solar values within the errors. The high- 
est oxygen abundance derived is 12-flog(0/II) = 8.85, 
about 1.5 solar0. The lowest oxygen abundance derived 



^The solar oxygen ab undance is taken as the value derived by 
lAspIund et aLl l|2005l ). 12+log(O/H)0 = 8.66±0.05 
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Fig. 6 Empirical calibration of the [SIII] electron temper- 
ature as a function of the abundance parameter SO23. The 
solid line represents a quadratic fit to the high metallicity 
HII region data 



is about 0.6 times solar. In all the observed CNSFR 
the 0/H abundance is dominated by the 0+/H+ con- 
tribution, as is also the case for high metallicity disc 
HII regions. For our observed regions, however also 
the S^/S^"*" ratio is larger than one. This is not the 
case for the high metallicity disc HII regions for which, 
in general, the sulphur abundances are dominated by 
S^+/H+. The derived N/0 ratios are in average larger 
than those found in high metallicity disc HII regions 
and they do not seem to follow the trend of N/0 vs 
O /H which marks the secondary behaviour of nitrogen. 
On the other hand, the S/0 ratios span a very narrow 
range between 0.6 and 0.8 of the solar value. 



5 Ionized gas characteristics 

From the calculated values of the number of Lyman a 
photons, Q(H'^), ionisation parameter and electron den- 
sity, it is possible to derive the size of the emitting re- 



gions as well as the filling factor (see lCastellanos et al 



20021 ). The derived sizes are between 1.5 arcsec and 5.7 
arcsec which correspond to linear dimensions between 
74 and 234 pc. The derived filling factors, between 6 x 
10~^ and 1 X 10~^, are lower than commonly found in 
giant HII regions (~ 0.01). Ha luminosities are larger 
than the typical ones found for disc HII regions and 
overlap with those measured in HII galaxies. The region 
with the largest Ha luminosity is R3-I-R4 in NGC 3504, 
for which a value of 2.02 x 10^" is measured. Ionising 
cluster masses range between 1.1 x 10^ and 4.7 x 10^ 
Mq but could be lower by factors between 1.5 and 15 
if the contribution by the underlying stellar population 
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Fig. 7 Distribution of derived oxygen abundances of 
CNSFR and high metallicity disc HII regions (top) as com- 
pared to that of the O23 parameter (bottom). While the 
first show comparable distributions for the two kinds of ob- 
jects, the second shows lower O23 (higher abundances) for 
CNSFR 



is taken into account. These values are consistent with 
those found from the kinematical data set. 

When compared to high metallicity disc HII regions, 
CNSFR show values of the O23 and the N2 parameters 
whose distributions are shifted to lower and higher val- 
ues respectively, hence, even though their derived oxy- 
gen and sulphur abundances are similar, higher values 
would in principle be obtained for the CNSFR if pure 
empirical methods were used to estimate abundances. 
This can be seen in Figure [7] where the distributions of 
the derived 0/H abundances and the O23 pararameter 
for CNSFR and high metallicity disc HII regions are 
shown. 

CNSFR also show lower ionisation parameters than 
their disc counterparts, as derived from the [SII]/[SIII] 
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Fig. 8 Distribution of tlie [SII]/[SIII] ratio for tlie observed 
CNSFR (dark) and tlie sample of high metallicity disc HII 
regions (light) showing the segregation between the two fam- 
ilies of objects in ionisation parameter 



ratio. Their ionisation structure also seems to be dif- 
ferent with CNSFR showing radiation field properties 
more similar to HII galaxies than to disc high metallic- 
ity HII regions. This can be seen in Figure[S]that shows 
[OII]/[OIII] ratio vs [SII]/[SIII]. Since the parameter r?, 
which is the quotient of these two ratios, is a go o d ind i- 
cator of ionising temperature IVilchez fc Pa"^ (|l988h . 
diagonal lines in this plot should run along objects with 
similar temperatures. CNSFR in this diagram also seg- 
regate from the high metallicity disc HII regions point- 
ing to higher values of the ionising temperature. 

The possible contamination of their spectra from 
hidden low luminosity AGN and/or shocks, as well as 
the probable presence of more than one velocity compo- 
nent in the ionised gas corresponding to kinematically 
distinct systems (see Section 3 above), should be fur- 
ther investigated. 
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